Communications with decode- and-forward relays in mesh networks 



Angel Bravo-Santos 
Universidad Carlos III de Madrid, 
abravo@tsc . ucSm . es 



Petar M. Djuric 
Stony Brook University 
djuricOece . sunysb . edu 



O 
(N 
+-> 
O 

O 



October 17, 2012 



Abstract 

We consider mesh networks composed of groups of relaying nodes which operate in decode- 
and-forward mode, where each node from a group relays information to all the nodes in the 
next group. We study these networks in two setups, one where the nodes have complete channel 
state information from the nodes that transmit to them, and another when they only have the 
statistics of the channel. We derive recursive expressions for the probabilities of errors of the 
nodes and present several implementations of detectors used in these networks. We compare the 
mesh networks with multihop networks, the latter being formed by a set of parallel sections of 
multiple relaying nodes. We demonstrate with numerous simulations that there are significant 
improvements in performance of mesh over multihop networks in various scenarios. 



> 

cn 
(N 

d 
X 



1 Introduction 



An old concept in radio communications, known as relaying, has been in the center of interest of 



various studies in communications Ephremides , 2002 , Goldsmith and Wicker 2002 , Ribeiro et al. 
[2008] . To a good extent, this interest has been driven by commercial applications and by the 
ability of wireless networks to exploit relaying so that they have reduced energy consumption and 
thereby increased lifetime. In general, relaying is used to provide for improved error performance 



and capacity Liu et al. , 2009 



In networks with relaying, nodes cooperate in moving information from source to destination iLiu 



et al. , 2009 . In the study of these networks, the important notion of cooperative diversity was 



introduced in A. Sendonaris et al. , 1998| , and then extended in A. Sendonaris et al. , 2003a and [j] 



N. Laneman et al. , 2004 . A class of wireless networks that are formed by a set of parallel sections 



of multiple relaying nodes was thoroughly studied in the literature Lin et al. , 2006 T. Wang et al. 



2007 . They are known as multi-branch multi-hop networks, and they can reach a full diversity 
order |T. Wang et al. 2007 . However, unlike mesh networks, these networks do not allow for full 
connectivity among nodes. 

Optimum detection in a network demands knowledge of all its channels. There is some work on this 
subject in the literature, and especially on particular network topologies. For example, in a recent 
paper, the problem of detection by a set of sensors and the communication of their information were 
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considered together [B. Chen et al. 2004 . Various fusion rules were proposed that correspond to 



different scenarios in terms of what is known for deriving the rules. The study was performed for 
multi-branch multi-hop topologies, but not for general mesh networks. 

An analysis of the maximum likelihood (ML) detector in cooperative multibranch multi-hop net- 
works was presented in ID. Chen and J. N. Laneman 2006 . A closed form expression was obtained 



for the bit error rate of a network with one relay. The difficulties with ML detection in DF net- 
works have led to several suboptimal alternatives. In A. Sendonaris et al. , 2003b , a variant of 



the maximal-ratio combiner (MRC), the A-MRC, was introduced to combine signals from several 
branches. In order to explicitly obtain a full diversity order, the cooperation MRC (CMRC) was 
proposed in T. Wang et al. , 2007 . For this detector, the weight of the channel of relay-destination 



was chosen to maximize the equivalent SNR for the channel. In that paper the multibranch multi- 
hop topology was considered. It is interesting to mention that a multi-node cooperative network 
can be viewed virtual MIMO system [Liu et aTj 



2009, A. K. Sadek et al., 2007 . In successive 



phases the relays combine the received signals from previous relays and the source using the MRC 
criterion. The system is not optimal and it can be adapted to a mesh network. 

Most of the above solutions assume perfect knowledge of the channel state information (CSI). 
However, having such knowledge can be expensive. For instance, in a sensor network the battery 
consumption can be high if the nodes have to inform about their CSI in a variable environment. In 
other circumstances it could be impossible to have complete knowledge of the CSI. For these cases 
it is important to develop schemes that operate with less information about the channel: instead of 
full information one has available only the channel statistics. There is some previous work in the 
literature on detection using channel statistics. In |A. Bravo et al. 2006 , the average SNR is used 



as a weight for combining signals in a multi-branch multi-hop relay network. In Lin et al. , 2006 



optimum ML fusion rules are derived for a joint sensor-communication problem with a multi-branch 
multi-hop topology. 

Data traffic in a sensor network is usually small. We take this statement as true, although it is not 
always the case, and we assume that most of the available channels can be used for cooperation and 



we do not put a limit to this number. This is a usual assumption in the literature, see T. Wang et al. 



2007 . A consequence of this assumption is that the spectral efficiency is reduced. In the literature, 
there are examples of simple networks where the error probability and the spectral efficiency are 
improved at the same time A. Sendonaris et al. , 2003b . However, it is not straightforward to 



extend these results to arbitrary networks. In this paper, we look for the objective of minimizing 
the transmitted power, or, what is related to it, minimizing the error probability, with no restrictions 
in the number of available channels. 

The cooperative detection problems mentioned above were posed for multibranch and/or multi-hop 
networks. For them, the optimal ML solution was abandoned in favor of tractable solutions. In this 
paper we present and analyze the optimal maximum a posteriori (MAP) detector for cooperative 
mesh networks, and we provide explicit analytical solution for a general topology. We consider 
optimal detection in two cases related to the channels in the network. One is when the CSI is 
available and the other when the channel statistics are only known. The first case was studied 
in part by the authors in [A. Bravo and P. M. Djuric 



20091. The relaying in the network is DF 



with uncoded and symbol-by-symbol demodulators. In summary, the two main contributions in 
our paper are the following: (a) we derive optimal MAP detectors for the nodes in mesh networks 



2 



and (b) we propose several implementations of the detectors for both, the case of known CSI and 
scenarios of known channel statistics. In the paper we do not consider the problems of routing and 
protocols that may arise, for example, in ad hoc mesh networks. How they relate to the proposed 
schemes here will be addressed elsewhere. 

The paper is organized as follows. In Section [2] we formulate the problem. In Section [3j we provide 
the general MAP solution that holds for known CSI as well as for known channel statistics. In 
Sections |4] and [5] we present detectors in the first and higher groups of nodes, respectively, and 
provide the specific solutions for the above two scenarios. Simulation results that demonstrate the 
performance of the mesh networks and how they compare with the multibranch and multi-hop 
networks are shown in Section [6j In the last Section [7j we have some concluding remarks about our 
findings. 



2 Problem statement 

We observe a mesh network where its nodes are grouped into relay groups and where the notation 

(k) (0) 

R- signifies the ith node of the fcth group. We denote the source by R\ and the destination by 
^(^+1)^ the superscript K + 1 implying that a message is relayed K times on its way from the source 

(3) 

to the destination. For example, R\ refers to the fourth node from the third group of nodes. We 
use the symbol n^^^ to denote the total number of relay nodes in group k. In Fig. [I] we show a 
drawing of a general mesh network of our interest. 




Figure 1: A general mesh network with K groups of relaying nodes. 

In this paper, for comparison purposes, we also work with multi-branch multihop networks. Multi- 
branch multihop networks can be obtained from mesh networks when we remove connections be- 
tween nodes so that each node in a group is connected to only one node from the previous group, 
as in Fig. [2| Clearly, a wireless network formed by a set of interconnected nodes can operate as 
a mesh network or as a multi-branch multihop network. In mesh networks, a node receives and 



3 



processes information from more than one node from a previous relay group, whereas in a multihop 
network, it does only from one node. 




Figure 2: A general multi-branch multi-hop network with K groups of relaying nodes. 
We consider binary modulations only (the extension to other modulations is straightforward), and 

(k) 

we assume phase coherent reception. The received signal by node and transmitted by node 
Rf' ^\ where A; > 1, is denoted by 

(k) Ak) (k-i) , (k) 

where G 1} is the transmitted symbol, h^j^ is the channel fading between i?-'^ and 

Rj^K and w^j^ is the channel noise. In this formulation, h\j^ > is the real- valued fading envelope 

of the channel. We denote the transmitted symbol by the source with x and not x^^^ so that we 
avoid superfluous notation. The noise is modeled as a zero mean Gaussian random variable with 
variance a^, i.e., w^j^ ~ A/'(0, cr^), and it is considered identically distributed in all the channels. 

We assume a simple protocol where every node has access to an orthogonal channel. This can be 
implemented with time division, frequency division or hybrid multiplexing. This assumption could 
be considered unrealistic, but with a little traffic offered in a network, as is frequently the case in 
sensor networks, most of the available channels can be used for cooperation with the node source 
of information and the time/frequency channels can be reused by distant nodes. We also point out 
that in a real system the number of available channels is not necessarily small. For example, the 



Zigbee standard for sensor networks ZigBee Alliance 2009 envisages a hybrid access for its next 



release in order to increase the number of channels and, thus, reduce message collisions. 
We distinguish two scenarios: 

1. all the channels in the network are completely known (CSI is available to all the nodes), and 

2. only the statistics of the channels are known, where it is assumed that the channels follow the 
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Rayleigh distribution, i.e., for the fading of a channel h we write 



2h - 

h - -^e ^h, h>0. (2) 



The signal-to-noise ratio (SNR) in the former case is defined by 



and in the latter, by 

„2 



7 = 32' (3) 



7 = ^- (4) 



Obviously, we can have a combination of the two scenarios, but we will not consider it here. 

Given the available information, the objective is to determine the optimal decisions at each node 
and study the performance of the systems under the two scenarios. 



3 General solution 

In this section we derive the decision rule of the nodes under the assumptions from the previous 

(k) 

section. We denote the probability of correct decision of node i in the kth group by ^ = 

P{xf^ = x\x) and the probability of error by pf" = P{xf = —x\x). In the sequel, we assume 

P(x = 1) = P{x = -1). 

Clearly, the nodes , j = 1, 2, • • • , n*^^) make their decisions according to 

' "p(x = -ii,<;')"/(!,i;v-=-i)' 

where j = 1,2, - ■ ■ , n^^\ and f{y^j\x = s) is the likelihood of x = s, with s being 1 or —1. 

The following groups of nodes will receive in general more than one signal and a decision is made 

(k) (k) 

by using all of them. Let the signals received by node Rj be denoted by the vector y^ = 
[yy\y2jV"" jyl^(k-i)jV 1 ^iid Ist the decisions of the nodes in the previous group, be given by 
x(fc-i) = [x'^^-'\x^t'\ ■ ■ ■ Then we have 



r(fc) 



^(x = l|yf) _ /(yf|x = l) 
P(x = -l|yf) /(yf |x = -l) 

Ex(^-i) /(yf |x('^-^))^(x('=-i)|x = -1) 



(6) 
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The decision rule at node Rj is 

xf = sgn(log(Lf )) . (7) 

For k > 2 and general mesh networks, the decision rule ([T]) requires that a node must know the 
states or statistics of the channels and the joint probability mass function (pmf) of correct decision 
of its relaying nodes, P(x(^~^)|x = 1) and P(x(^~^)|x = —1). The former requirement in most 
studies is assumed satisfied. The second requirement is much more restrictive, and in the sequel we 
propose ways of dealing with it. 

We note that in ([g]), we have likelihoods of the form /(y^^^^lx). For them we have the following 
claim: 

Claim 1: For the likelihood /(yj'^''|a;), A; = 1, 2 in mesh networks, we can write 

/(yfk) = 11/ (4^1-)- (8) 

1=1 



Proof: The proof follows directly from the conditional independence of the channels and the decision 
independence of the nodes in the first group. 



From Claim 1, for the likelihood ratio in nodes of group 1 and 2, L^j\ k = 1,2; j = 1,2, ■ ■ ■ , n^'^)- 
we can write 



if = n if = n !^ 



where 



= = 1\X = 1) = P(xf-') = -1\X = -1) (10) 

is the probability of correct decision, and 

pf = 1 - (11) 
is the probability of error, with P(0) = 1. 

The likelihood ratio in ([9]) , in general, does not have the same simple form for k > 2 because of the 
dependence among the decisions of the relaying nodes. However, if we approximate 

„{k-i) 



P(x('=-i)|x) « H P(xt'^\x), (12) 



i=l 



we will obtain suboptimal detectors that use ([9]), as is shown in Section [5j We note that with the 
assumption ( 12 ), we only need to know the marginal probabilities of correct decisions of the relaying 
nodes. 
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A symmetry property, similar to (10) and (11), can be stated for P(x^'^ ^)|a; = 1) and P(yS-^ 



X 



— 1). We make the following claim: 

C/aim ^: If the channel likelihoods satisfy /(y^^^"* I x-'^ = /(~2/ij^l ^'^'^ f {,~y^if\^^t ^'') ~ 

/(^ij^l ~ ^o'^ ^' have 

P(x('=)|x = l) =p(-xW|x = -l). (13) 

Proof: See Appendix 1. 

We note that if the conditions for Claim 2 are satisfied, for optimal processing of the received signals 
all the information about the joint pmf of correct decisions needed by a node in the kih. group is in 
p^^(fe-i)|^ _ Then, for the likelihood ratio ([6| we can formally write 



Ex(-i)^(-x('=-i)k = i)n£r'/(ylfix.(^'-i))' 



where all the probabilities are conditioned on x = 1. 



4 Detectors in the first group of nodes 

In Section [2j we described two scenarios, one where the necessary CSI is available to the nodes, 
and another where only statistics of the channels are known. Here we describe the detectors of the 
nodes in the first group and we present their performances. 



4.1 Completely known channels 

In the case of known channels and based on the assumptions from Section [2| we can express the 
likelihood as follows: 



The decision rule of node R^P is based on 



p{yij\x = -1) 



p{yij\x = 1) 



which simplifies to 



xf = sgnfyjjM. (17) 



For the nodes in this group, we can easily find the probability of error, and it is given by 
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with Qiz) = 1 — where ^{z) = ~^2"^ ^ ^^"^ ^ij^ SNR of the channel that hnks 

the transmitter and the j'th node in group one (see ([3| 

4.2 Channels with known statistics 

When the channels are not known, and instead the nodes only have available the channel statistics, 
we proceed as follows [Niu et ah 2006 . We write for the likelihoods 



= / f{^\f^.^fih^)d^, (19) 



where f{hfj) is given by Q. 



For simplicity, we rewrite the last integral without subscripts and superscripts, and we get 
It is easy to show that f{y\x) can further be expressed as 



f[y\x) = ^-J-^ exp(^-^^ ^ hexp{-fih^ + 2uh)dh (20) 

where 

1 y X 

^=2^ + ^'"=2^- ^'^^ 
We can show that the integral in the above equation can be analytically solved. We get 



After substituting the solution (22) in (20) and returning all the superscripts and subscripts, we 
obtain 

where 




In Fig. [3| we plotted the probability density functions (pdfs) of y given x = ±1 for SNR = 3dB 
for the two cases of interest, when the channels are completely known and when their statistics are 
only known. As expected, the latter case leads to deterioration in performance because then the 
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0.7 




Figure 3: The pdfs of y, f{y\x), for x = ±1 when the channels are known (curves with dash-dot 
hues) and when the statistics of the channels are known (solid lines) for SNR equal to 3dB. 

pdfs are shifted towards y = and are flatter in comparison to the pdfs that correspond to known 
channels. 

For the probability of error at node Rj^^ , we have 

pf^ = p (l'i^ < 1|x = l) = P{y <0\x = l), (24) 



where 



I.e., 



= /°^/f^exp(-^)(l + V^a,$(a,)exp(^^))d,. (25) 
Note th at P^^^ i s the probability of error of a coherent binary signaling over a Rayleigh fading 



channel [Proakis 2001 . The above integral can be solved and expressed in a closed form by using 
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integration by parts, and the result is given by 




5 Detectors in the groups of nodes beyond the first group 



The detectors for the nodes in groups k > 1 are challenging to implement because of the need to 
know the joint probabilities of correct decisions of the relaying nodes (see Q). It is interesting to 
point out that there is one exception to the need for pmfs, and it is the case of mesh networks with 
two nodes per group, known CSI, and assumptions as in Claim 2. Here, we describe four types of 
detectors for the nodes Ri'\ A; > 1, i = 1, 2, • • • ,n^^\ We also briefly describe the detectors in 
multihop networks. 



5.1 Estimation of the joint pmf based on Monte Carlo sampling (MCS) 

The joint pmf of interest is P (x'^'^-'lx), and we can show that it can be expressed as 

P (x('=)|x) = P (sgn (logLS'^)) , sgn (logL^^^^) , • • • , sgn (logL^';!)) |x) , (27) 

(k) 

where the random variables are not independent. 

From ([6]), we see that the L^j'^s are expressed via the decisions a;'^'^"^^ which allows us to propose 
a computational method for obtaining the joint pmf. The procedure is the following: 

1. Each node ijj'^'' generates samples where i = 1,2,- •• ,n^^~^\ m = 1,2, ••• ,M from 
the mixture distribution given by 

„{fc-i) 

/(yi„?/2,v,y,.v,y„(.-i),) = E ^(-^'^'^1^) n 

(k) (h\ (kY^^ 

2. The node computes the signs of the obtained logL^- , i = 1, 2, • • • , ' at the drawn ylj , 
J = 1, 2, • • • , n^^^ and stores them as a vector in a counter. 

3. The process is repeated M times. 

4. From the obtained outcomes, the node obtains the estimates of the joint probabilities. 



The generation of samples in the case of a mixture Gaussian is easy. When the channel statistics 



are known only, the sampling of yij must come from a distribution given by (23). To that end, one 
could apply rejection sampling. 

Clearly, at the end, all the nodes in the kth group will have different estimates of the joint pmf. The 
next group may combine all these estimates, for example, by taking the average of the estimates. 
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This approach becomes tedious when the number of nodes in a group gets large. For example, when 
n^*^) = 10, one has to work with a joint pmf of a 10-dimensional vector, which, in general, requires 
the computation of 1023 estimates. However, there are two ways to decrease the computational 
load of this approach and they are (1) to have a designated node in a group compute the joint pmf 
and (2) under an assumption of symmetry, instead of computing the values of 2"- — 1 elements, a 
node computes only n*^'^-* + 1 elements. 

5.2 Estimation of the joint pmf based on pilot signals (PS) 

Suppose that node Rf^K k > 2 receives messages from n*^*^"^) relay nodes, y^j\ In the rest of the 
subsection, we suppress the subscript j and the superscripts (k). For example, for the received 

(k) 

signals ylj , we write y^. 

Define the variable = 1 if j/j > and = if < 0. Let also 

n n 

n{zi,Z2,--- ,zn) = ^Zi2^-^ = Xl'^i2^-\ (29) 

i=l i=l 

where Ki € {0, 1} identifies the ith binary digit of k represented in binary notation. Thus, for 
example, k = corresponds to {yi < 0, j/2 < 0, • • • , y„ < 0} and k = 1 to {yi > 0, ?/2 < 0, • • • , j/„ < 
0}. In other words, the joint pmf of zi,Z2, ■ ■ ■ ,Zn can succinctly be represented by the pmf of k . 

Similarly, we represent the joint pmf of the transmitted symbols of the relaying nodes. Let ( and 
(i be defined by 

n n 

C{xi,X2,--- ,Xn) = J2{xi + l)2'-^ = J2Ci^'~'^ ae{0,l}, (30) 

i=l i=l 

and the joint pmf of Xi is equivalent to the pmf of (. 

We want to estimate p{C) from pilot data (where the source, say, repeatedly transmits x = 1). Then 
the node first estimates the pmf of k and, based on the model for the channels, estimates the pmf 
of (. We can write for «; = 0, 1, • • • , 2" — 1 

2"-l 

p{k) = J^p(k|c = OKC = 0, (31) 

1=0 

where the probabilities p{k\(^ = I) are known, i.e., we can show that 

n 

Pi^\c=i) = n^'cf Hi'^^'^''^ (32) 

i=l 

where 
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where Ki = Q means 

_ Xi + 1 

and Pc^i is the probabihty of correct decision based on the transmitted signal from the ith node, 
and Pe,i is the corresponding probabihty of error, that is pe,i = 1 — Pc,i- This probabihty of error is 
obtained from (26). 

We can represent the probabihties of k by 

p« = [p{k = 0)p{k = 1) ■■■ p{k = 2" - 1)]^ 
and the probabihties of C by 

PC = [p(C = 0) p(C = 1) ••• p{C = 2^-i)V. 



We can also construct a 2" x 2" matrix, P with probabilities as the ones defined by ( 32 ) . Then we 
can write 

= PPC- (34) 

We note that the matrix P is a doubly stochastic matrix, that is, its rows and columns sum up to 
one. 

We reiterate that we estimate the vector from the pilot data, and we denote it by p^- Then the 
estimate of p^ is readily obtained from 

PC = P^'Pk- (35) 

This estimate, however, may produce negative probabilities (although the sum of all the elements 
of P(^ equals one), which is due to errors in estimating p^. Thus, we change the set of equations 
(34) to the following minimization problem with linear constraints: 

p, = P(pC + e), (36) 

where e has minimum norm and adds to 0. In Appendix 2, we show that the solution vector p^ to 
the system ( 36 ) has elements given by 

[Pc]^ = ib^-0^, (37) 

where (x)"*" is the positive part of x, z = 1, 2, • • • , 2", b = P^^p^, and ^ is a parameter related to 
the Lagrange multipliers. It is worth pointing out that (37) is in fact an upside-down waterfilling 
scheme. 

In summary, according to this scheme the source node transmits symbols x known to the rest of 
the network. The nodes in groups k > 1 can estimate from the received signals from their relaying 
groups the joint pmfs P(x('^~"'^)|x) as just shown. When the estimation is completed, the nodes 
have the necessary information for processing of the signals with unknown transmitted symbols. 
Obviously, a drawback of the scheme is that some transmitted symbols are used for estimating joint 
pmfs and some power is used for estimating the pmfs. 
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5.3 Estimation of the joint pmf based on a predefined joint pmf (PJP) 



The biggest values of the joint pmf Pyx\ 



(fe-i) Jk-i) 



= l) are for combinations when 



most of the arguments x 



(fe-i) 



1. This suggests that we approximate the joint pmf by 



P(x 



(k-i) Jk-i) 



(fc-i) I 



1) 



\{i:x\ 
otherwise 



l,iG{l,2,. 



(fc-i) 



}}\>Nf 



, (38) 



where pc is probabihty obtained from the number of combinations of the elements of x'-'^ that 
satisfy the condition in (38), |{-}| is the cardinality of the set, and Nf < n^^~^\ In other words, the 



probability of correct decision is constant and different from zero whenever the number of arguments 
^{k 1) _ ij^g enough. This detector is efficient because only a few terms in the numerator of 
the likelihood ratio (6) are different from 0. Claim 2 is used for selecting the significant terms of 
the denominator of (6). This simple detector does not require the use of pilot signals. 



5.4 Estimation of the joint pmf based on the assumption of independent deci- 
sions (ID) 



When for k > 2 we approximate the joint pmf of x^^ by ( 12 ), for decision making we can use the 
hkehhood ratio (|9|. If we take the logarithm on both sides of ([9]), and we define 

= log LSf = log (L^^Pt'' + Pf') - log (Pf' + Pf') , (39) 



where 



l^f = logf{yi;'\x = l)-logf{yi;'\x = -l), (40) 



and 



" /(i/rkr^) = -i) ^ ^ 

(k—l) 

representing the likelihood ratio for the transmitted symbol of i?] , we can write for the overall 
loglikelihood ratio 

1^ = Y.^- (42) 



i=l 



The decision rule then simplifies to 



xf> = sgn fM. (43) 
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(k) 

For the probability of error at Rj , we can write 



pf = p /W<o|x=l =P j;/W<0|x = l 



j=i 



log ( + pf^) - log (f^'^-^) + 



< X = 1 



(44) 



(k) 

We note that under these assumptions the new random variables l^j are independent. Clearly, if 

(k) (k) 

the Uij s are identically distributed, so are the loglikelihood ratios 1- . 

The implementation of this detector requires information about the probability of correct decision 
of the nodes in the previous group P-'' ^\ i = I,-- - ,n^'^^^'^. These probabilities can readily be 
computed as described in subsections A and B. The computations in this case are much simpler 
because the problem is broken into n^^~^') independent problems, where in each problem we estimate 
the probability of detection (error) of only one relaying node at a time. 



Another possibility is that the nodes compute their own probabilities of errors from (44) and trans- 
mit them to the nodes of the next group. The advantage of this approach is that there is no need 
for transmitting pilot signals. We note that in the transmission of the probability, some form of 
quantization must be adopted. In Section [6j we study the influence of quantization via simulations. 

We now briefly describe the implementation of the scheme in the cases of completely known channels 
and known channel statistics, respectively. 

5.4.1 Completely knovi^n channels 

For the conditional distributions of the observations, we can write 



and for the likelihood ratios. 



/(l/fk = l) = Pt'^M{ff,a^)+pt'^M{-f^\a^) 



(45) 



4« = expf?^). (46) 



whereas for the loglikelihood terms we have 



4f = logfexpf?^U(-^) + pK 



log I Pf^ + exp ( ^1^=4^ 1 Pr^' 1 . (47) 



a2 
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(k) 

We can also find the distribution of /)^- by change of variables. We obtain 



p{k-i)^ '^.i . (48) 



1/ . 1 pC^-D-pC^-De^'f hlj\2 



p(^-l)^ ^J P^'^-De'-J -P.C^-l) 

+ -T,- e s » 



The probabilities of error of the nodes in group k > 2 are obtained by the recursive equation (44). 
A closed form analytical solution of the recursive equation in the general case, however, cannot be 
obtained. 



5.4.2 Channels with knovi^n statistics 



The probability of error at the nodes in group A; > 2 is the same, that is, it is given by (44), where 
now 



L 



(k) 



l-V2^ayfcI>(-ay(f)expf^^ 



(49) 



is obtained by using (41 ) and (23). The probabilities of error of the nodes in groups k > 2 are again 



obtained numerically by the recursive equation (44). 



5.5 Detectors for multihop networks 



In multihop networks, the variables y. 



(fc) 



1,2,- 



, n 



k > 1, are conditionally independent. 



Therefore, for detection at the relay nodes one uses the likelihood ratio in The last node in 



the network, i?[ , collects the information of all the nodes from the previous group and uses 

. — . (f,\ . — . 

the decision rule (43), where the sum loglikelihood Ij is calculated from (42). For this node, the 

obtained results under the independence assumption apply. 



6 Simulations 



We conducted many experiments where we compared the performances of the mesh and multihop 
networks shown in Figs. [T] and [2| respectively. In the experiments, the channels were Rayleigh 
distributed and, for better understanding of the results and easier comparisons, the SNRs were the 
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Figure 4: Probabilities of error for mesh and multihop networks as functions of number hops for 
the case of known statistics. 



same for all the nodes. Most of the simulations were for the scenario where the channel statistics 
are known. This case presents two main advantages in real networks: (a) the information about 
the statistics is provided just once for each node and (b), the probabilities of correct detection are 
evaluated also once. By contrast, when detecting with complete CSI at the nodes, the channel 
information must be updated if the CSI changes, and the probability of correct detection has to be 
evaluated accordingly. 

First, we considered a mesh network with a decision rule based on complete knowledge of the 
probabilities of correct decision of the nodes of the previous groups and based on ([T]) . The results 
of this setup were a benchmark for comparison. The joint pmfs of correct /incorrect decisions of the 
groups for k > 1 were computed using the MCS method. 

In Fig. [4j we show in the curve Mesh MCS the probabilities of error of mesh networks with 
different number of hops when all the joint pmfs were known. In the simulations, the number of 
nodes per group was 10 and the SNR = 3 dB. The error probability decreased as the number of 
hops increased until a point where it remained constant with the number of hops. The spatial 
redundancy introduced by the mesh network lowered the error probability until an error floor was 
reached. 
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In Fig. [4j we also present the performance of the mesh network using the PS method. As can be 
seen, the obtained results are close to the ones of mesh networks that have complete information 
about the pmfs. 

To these curves, four other curves are also displayed. One of them shows the performance of the P JP 
detector. With the curve Mesh ID Q=4 bits, we show the performance of the detector based on the 
ID detector, and where the probabilities of correct decision are quantized with four bits. Note that 
a node uses the probabilities of correct decision of previous nodes and quantizing these variables 
reduces the amount of information to be distributed. Finally, with the Multihop ID Q=4 bits curve 
we plotted the performance of the multihop network, and with the MRC curve, we displayed the 
performance of the MRC detector [Proakis 2001 . For the multihop network, the probabilities of 



correct decision of the previous node were again quantized with four bits. The MRC detector was 
implemented by combining the signals arriving to a node using the MRC criterion. 

The worst performance was clearly achieved by the multihop network. Also, its performance dete- 
riorated steadily as the number of hops increased. The performance of the network whose nodes 
had to estimate the joint decision pmfs of the nodes in the previous group performed almost as well 
as the network whose nodes knew the joint pmfs. Surprisingly well performed the mesh network 
that used the independence assumption and where the nodes quantized their probabilities with four 
bits. Next in performance of the mesh networks came the one that employed the simple detector 
and finally, the worst was the network with MRC detectors. 

The analytical study of the quantization effect on the ID detectors is difficult. Instead, we performed 
experiments where the number of bits for representing the quantized probabilities was either one 
bit (the quantized probabilities were zero and one) or four bits (for 16 probability values, that is, 
the quantized probabilities were 0, 0.0666, • • • , 1). We also obtained results when the probabilities 
were not quantized. The results are shown in Fig. [5] . There we see the probability of error as a 
function of the number of hops. The nodes assumed independence of the decisions of the nodes of 
the previous group. The SNR was 3 dB. In all the mesh networks with 10 nodes per group, the 
BER decreased with the number of hops. The results show that four bits were enough for achieving 
almost as good performance as the one when there was no quantization. Interestingly, using only 
one bit for quantization also yielded good results. Note that representing the probability of correct 
decision with one bit is a form of selection combining: the nodes with low probability withdraw 
themselves from combining in the following group of nodes. 

In the case of five nodes per group (curve Mesh ID 5 nods/grp Q=4 bits) the decrease in probability 
of error probability was slower. The performance of this network was worse than that with 10 nodes 
per group by an order of magnitude. In Fig. [5| we also plotted the performance of two multihop 
networks. One of them was the same network whose performance is shown in Fig. |4j and the other 
corresponds to the multihop network with CMRC detectors (Multihop CMRC 10 nods/grp) fT^ 



Wang et al. , 2007 . Note that the CMRC detector uses complete information about the channel. 
Nevertheless, it has a larger error probability than the multihop ID detector (curve Multihop ID 
10 nods/grp Q=4 bits), where only channel statistics are used. In both cases the error probability 
increases very fast with the number of hops. 

Finally, in Fig. [6] we present curves of probability of error versus SNR for a mesh network with 
MCS detectors, a mesh network with ID detectors and four bits for quantization, and a multihop 
network. The number of nodes per group was 10 and the curves were obtained after 9 hops. 
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Figure 5: Probabilities of error for mesh and multihop networks as functions of number hops for 
MAP detectors that used conditional independence among the probabilities of correct decision. The 
nodes knew the statistics of the channels and the SNR was 3dB. 



Prom the experiment we can draw the following conclusions: a) the assumption of independence is 
adequate and gives results that are very close to the case of dependence between the probabilities 
of correct decision of the nodes in a given group and b) the mesh networks offer great advantages 
in performance over multihop networks. 



7 Conclusions 

In this paper we studied the performance of mesh networks whose nodes operate in the decode- 
and-forward mode. In particular, we investigated the probability of error in these networks as a 
function of SNR, number of nodes per group, and number of hops. Wc studied two cases, where in 
the first case the nodes knew the CSI and in the second, the nodes worked only with the statistics 
of the channels. We compared the mesh networks with multihop networks and showed the gain 
in performance of the former with respect to the latter. In the presented work we dealt with 
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SNRdB 

Figure 6: Probabilities of error for mesh and multihop networks as functions of SNR for MAP 
detectors. 

binary modulations, but the generalization of the proposed detectors to deal with M— ary, M > 2, 
modulations is not difficult. 



Appendix 1 
Proof of Claim 2 

We prove the claim by induction. Suppose that the channel likelihoods satisfy f{ylj^\x['' = 
/( - Vij^l - ~^^) and /( - yij''\xf~^^) = fiVij^l - x^'''^''), and that the claim is true for k - 1, 
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i,e. P(x(''-i)|x = 1) = P(-x('=-^)|x = -1) Then we have the hkehhood ratio (|6j) 



' Ex(-^)nr=r/(yifix.('=-i))p(x(^-i)k = -i) 
Ex(^-^) nrrr^' /(-4V:r.(^-i))p(-x(^-i)|. = -i) 
Ex(-i) nil"' /(-ySf i-x.(^-i))p(-x(^-i)|x = 1) 



where we have used the inductive hypothesis and the symmetry of the Ukehhoods. We can reorder 
the sums by replacing x with — x to get 

_ Ex(-^) nir ixf-^))p(x('=-i)|. = -1) 



1 



(51) 



From (51), we have the equahty 



^'(^r (-yD 1 1> • • • , ^rc-yf) 1 1, • • • |x = 1) . (52) 

Next we show that 

/(-yfk = i)= E n /(-ylf|x,('=-^))P(x('=-i)|x = i) 

x(fc-i) «=i 

= /(yf k = -i) (53) 

and therefore 

P(Lf)(yf))|l,...,Lf)(yf)|l,...|x = l) = 

P{L? iy['^) I 1, • • • , i^f (yf ) I 1, • • • k = -1) (54) 

proving the claim for the kth step. 

The proof of the claim for /c = 1 is immediate. This completes the proof of Claim 2. 



Appendix 2 



The system (36) can be expressed as the following optimization problem: 
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minimize 



E 



subject to pf + €i 



bi, i = 1, 



(fc-i) 



-pf-i) <0, i = l,.-- ,n(^-i) 



(55) 



The last two conditions make sure that the vector pC^"^) is a probabiHty vector because, as can 
easily be proved, b is a vector whose elements add to one. The problem (55) is convex and for it 
the Karush-Kuhn- Tucker (KKT) condition is necessary and sufficient for the optimal solution 



Luo and W. Yu, 2006 . More specifically, the KKT condition is 

„(fc-i) 



n(fc-i) 
i=l 



(fc-1) 



Pi 



< i = 1,- 

Pl ' +ei = bi i = 1, 

E^^ = o 

0, A,; > 



1)+ ^ r/iV(pf ''^ + e^-bi) + UeV^ei 

i=l 
(fe-1) 



n 



(k-i) 



(56) 



H-i^-'' - 1) 



; = i,...,n(^-i). 

After applying the partial derivatives, for the KKT condition we obtain 



+ ^(^i + ^e) = 0; 
Ai = J^j; 

Ai>0,A, = Oifpf~') >0, 

(fc-i) ^ 
K > 0, 



(fc-i) 



(fe-i) 
(fc-i) 
(fc-i) 



(57) 



(fc-i) 



Pi 



+ ei = bi 



(fc-i) 



We define the variable 



M 



(58) 



where M is the number of variables p^ equal to zero. With it and the KKT condition (57), the 
solution to the system ( 55 ) can be expressed in a compact way as given by ( 37 ) . 
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